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Abstract 

A deeper understanding of the vacuum structure in QCD invites one to rethink 
certain aspects of the axion physics. The recent advances are mostly due to devel- 
opments in supersymmetric gauge theories and the brane theory, in which QCD can 
be embedded. They include, but are not limited to, the studies of metastable vacua 
in multicolor gluodynamics, and the domain walls. We briefly review basics of the 
axion physics and then present a modern perspective on a rich interplay between 
the QCD vacuum structure and axion physics. 



An expanded version of a review talk given by G.G. at "Continuous Advances in 
QCD 2002/ ' Arkady] "est", honoring the 60'th birthday of Arkady Vainshtein, Min- 
neapolis, May 17 - May 23, 2002. 



1 Introduction 



Almost 25 years elapsed since the axion was introduced in particle physics |], ||] as a 
possible solution of the strong CP problem. Since then, it became a text-book and 
encyclopedia subject. For instance, Oxford University's "A Dictionary of Science" 
defines axion as "a hypothetical elementary particle postulated to explain why there 
is no observed CP violation (see CP invariance) in the strong interaction (see fun- 
damental interactions). Axions have not been detected experimentally, although it 
has been possible to put limits on their mass and other properties from the effects 
that they would have on some astrophysical phenomena (e.g. the cooling of stars). 
It has also been suggested that they may account for some or all of the missing mass 
in the universe." 

While at the early stages the axion physics was considered predominantly in 
the context of quantum chromodynamics, at present the center of gravity of the 
axion studies shifted in astrophysics. It was realized rather early that the axion was 
a viable dark matter candidate. The research on this aspect of the axion physics 
quickly picked up and never subsided. Extensive investigations were and are being 
carried out in the astrophysical community. At the same time, after the rapid 
advances in the 1980's, the QCD practitioners seemingly lost interest in this subject. 
The reason is obvious: the progress in understanding the QCD vacuum structure 
was painfully slow. The prevailing impression was that "nothing happened here," 
so there was no motivation for revisiting QCD-related aspects of the axion physics. 

In this review we will try to argue that "something interesting happened here." 
A substantial progress has been achieved in the recent years mainly due to insights 
in QCD obtained from supersymmetry and the brane theory. The existence of a 
multitude of (quasi) stable vacua at large N c and "abnormally" thin domain walls 
with "abnormally" heavy excitations are just a few topics of interest that should be 
mentioned in this context. A deeper understanding of the QCD vacuum structure 
requires a reassessment of a number of issues of direct relevance to axions. After a 
brief summary of basics of the axion physics we review these new developments. 



2 The strong CP problem 
2.1 The 6 term 

After the discovery of asymptotic freedom in QCD in 1973 [Q, Q for a short while 
it was believed that QCD possesses the same natural conservation laws as its more 
primitive predecessor, QED. The discovery that P and T conservation in QCD is 
not natural came as a shocking surprise. This fact was realized with the advent of 
instantons []5J which demonstrated that the so-called 9 term 

6 

ACe = ^^ G %G%, (1) 
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does not necessarily vanish. Here the dual field strength is defined as 

= ( V 2 ) £ nva/3 Gap • 

(The indices are assumed to be contracted via the flat space metric). The operator 
GG has dimension four, it can and should be added to the QCD Lagrangian. The 
9 term (jl]) violates P and T invariance (and hence, it violates CP since CPT is 
preserved). Note that the analogous term A £ ~ FF in QED has no impact on the 
theory whatsoever. What is the difference? 
The 9 term can be rewritten as follows 

A£ e = 9d ll K lt , (2) 

where is the Chern-Simons current defined as 

K/i = Yg~^2 £ V va P C^t 9a A a p + 2 fabcA-l A a A°p \ . (3) 

Being a total derivative, the 9 term does not affect the equations of motion. At a 
naive level, one can discard in the action the integrals over full derivatives. This 
was a rationale behind the original belief that QCD naturally conserves P and T. 

The instantons revealed the fact that the vacuum structure in QCD is more com- 
plicated than that in QED. In particular, the field configurations with the instanton 
boundary conditions give rise to a nonvanishing 

(ASg) onc inst = J d 4 X (A £,j)one inst = 9 . (4) 

The integral over the full derivative does not vanish []. Therefore, CP- violating 
effects may be present in strong interactions. 

We pause here to make an explanatory remark regarding Eq. (Q). A key notion 
is the topological charge V of a gauge field configuration, 

V = J d 4 xd,K, = J d'xK^x^WZt™ 

ee K(t = +00) - K(t = -00) , (5) 

where K, is usually referred to as the Pontryagin number. The topological charge is 
zero for any perturbative gauge fields — such fields are said to have trivial topology. 
The instanton field configuration has a nontrivial topology. In the Ao = gauge it 
interpolates between A m (x, t — > — oo) = , m =, 1, 2, 3, and 

A m (x,t -»• +oo) = U + 8 m U, (6) 

1 We jumped here from the Minkowski to the Euclidean formulation of the theory. In passing 
from Minkowski to Euclidean, the term (|l|) acquires an "i" factor, so does the integration measure 
in the action. Since this is a text-book topic, we will pass freely from Minkowski to Euclidean and 
back making no explicit statements as to which space any given formula belongs to. 
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where the matrix U is Poyakov's hedgehog, 



\/x l + p 

The Pontryagin number for U reduces to 



m = «p ( -^ra • ( 7 ) 



K. = — ^ | rf 3 x e iifc Tr ([/+ 9* [/) (U + & U) (U + d k u)=l, (8) 
implying that the instanton topological charge is unity. 

2.2 Superselection rule and 6 sectors 

As we pointed out in the previous section, the value of the parameter 9 a priori can 
be arbitrary. The theories with different values of 9 describe different worlds which 
do not "communicate" with each other. In other words, the worlds with different 9 
belong to distinct superselection sectors || [7J. 

To see that this is the case and to further elucidate the role of 9 let us consider 
pure gluodynamics in the Hamilton gauge A = 0. In this gauge the Lagrangian 
does not depend on A and the Gauss' law (which in other gauges could have been 
obtained by varying the action with respect to Aq) is imposed as a constraint on 
physical states 

D^iolPhys) = 0. (9) 

This gauge fixing does not eliminate, however, the gauge freedom completely. Purely 
spatial gauge transformations (independent of the time variable) are still allowed. 
The generator of these residual gauge transformations can be written as 

G(a) = exp(z j rf 3 xTr^G i0 (^)a(x)) , (10) 

where a = a a t a and the trace Tr runs over the color indices. The generator Q (a) 
acts on the spatial components of the gauge fields, 

g + A k g = u + (A k + id k ) u, (ii) 

where U = exp(ia(x)). Furthermore, it is straightforward to show that the operator 
Q does not commute with V if the corresponding gauge transformations give rise to a 
nonzero right-hand side in Eq. @ (the latter are called large gauge transformations), 

[QV] + 0. (12) 

Therefore, an eigenstate of V cannot be a physical state. Instead, the physical state 
is defined as a superposition of the eigenstates of V 

\9) = £ e * e », (13) 
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where V\n) = n\n). 

In other words, in the infinitely-dimensional space of fields there is one direction 
parametrized by the variable K, = J K d 3 x which forms a closed circle. The wave 
function (as a function of /C) is the Bloch superposition (|T3|) . The parameter 9 is 
nothing but a "quasimomentum" || BL f7|. In QCD it is called the vacuum angle. 

In this formulation the 9 angle enters as an arbitrary phase in ([13|) and varies in 
the interval from to 2tt. Physics must be 2tt periodic in 9. At 9 ^ or it one can 
expect P and T noninvariant effects. 

It is straightforward to show that for any gauge invariant operator O 

(9'\d\9) ~ cW (14) 

Therefore, no gauge invariant operator can transform a state of one 9 world into a 
state of another 9 world. The different 9 worlds are disconnected from each other. 

2.3 Constraints on 

As was mentioned, nonvanishing 9 leads to CP violating observables in QCD. (The 
point 9 = 7T will be discussed separately). It is known that strong interactions 
conserve CP. Hence, a natural question arises as to what are the experimental 
constraints on the value of 9. 

In a full theory, with quarks, there is an additional contribution to the CP odd 
part of the Lagrangian. It comes from the imaginary phases of the quark mass 
matrix A4. These phases can be rotated away from the mass matrix by chiral 
transformations of the quark fields. However, because of the axial anomaly ||, 
which manifests itself as a noninvariance of the Feynman integral measure under 
the chiral transformations at the quantum level [JTOj , the quark mass matrix phases 
appear in front of the GG term in the QCD Lagrangian. Therefore, the actual 
parameter that sets the magnitude of CP violation in QCD is 

9 + arg(det M) . (15) 

In what follows, for simplicity of notation, we will denote this parameter by the 
same letter 9, implying that the part arg(detA^) is included by default. 

Perhaps, the most pronounced effect of the 9 term is generation of a nonzero 
neutron electric dipole moment (nEDM). The latter is parametrized by the following 
effective Lagrangian 

(p _ 

CnEDM = ^175^™^' ( 16 ) 

where the photon field strength is F^ u = d^A u — d u A^, and n stands for the neu- 
tron. Moreover, is the antisymmetric product of two Dirac's gamma matrices, 
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In the presence of the 9 term, nEDM can be found from the following matrix 
element: 



(nip^^leJ^A^ ■ i J d 4 xA£ e \n( Pi )) = d? n n(jp } ) i 5 a^n( Pl ) k^(k) , (17) 

where J^ m is the quark electromagnetic current. The momentum carried by the 
photon, k^ = ipf)fi — {pi)n, equals to the difference between the final and initial 
momenta of the neutron while e u {k) denotes the photon polarization four- vector. 

The matrix element on the left-hand side of (|17| ) is a highly nonperturbative 
object. Its calculation in QCD is nontrivial. Nevertheless, there are a number of 
different methods by which nEDM had been estimated in the past. We will list 



them below. The bag model calculation was performed in Ref. |TTJ. The result is 



d^Bag - 2.7-l(r 16 e-cm. (18) 
Shortly after Ref. JTTI], the chiral logarithms (CL) method was used (T2|] leading to 



the estimate 

dl\cL - 5.2- l(T 16 e-cm. (19) 
The method of chiral perturbation theory (ChPT) was further advanced in Ref. (13 



with the following result 

G^|chPT ^ 9 3.3 • 1(T 16 e • cm . (20) 



Finally, the most recent paper on QCD sum rule (SR) calculations |14] of d^ n gives 



dl\sR ~ 9 1.2- l(T 16 e-cm. (21) 

All results above have a considerable uncertainty, of at least 50%, which reflects a 
variety of uncertainties inherent to nonperturbative QCD calculations. Even though 
the results scatter by a factor of several units, it is beyond any doubt that d1\ t h eor ~ 
#10- 16 e-cm. 

This number should be compared with the most recent experimental result for 
nEDM presented in Ref. |16| 

|^| exp < 6.3- 10" 26 e-cm. (22) 

One gets a very strong constraint on the value of 9, 

\9\ £KT 9 . (23) 

We see that if the value of 9 is nonvanishing it has to be unnaturally small. 
There is no a priori reason why two terms in fll5|), the bare theta parameter and 
arg (det Ai), should cancel each other with such an extraordinary accuracy, of one 
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part in 10 9 or better. A dynamical mechanism is needed to explain the unnatural 
smallness of the 9 term. 

Before proceeding further, let us mention that other CP odd effects are induced 
by the 9 term too. They impose less stringent bounds on 9, however. For instance, a 
nonzero 9 gives rise to a nonvanishing amplitude of a CP violating decay rj — > 7r + 7r~ 
for which Bt(i] — > 7r + 7r~) ~ 9 2 1.8 ■ 10 2 ]17j [U|. The experimental limit for this 
decay is Br exp (?7 — > 7r + 7r _ ) < 1.5 ■ 10~ 3 . This yields a constraint \9\ < 3 • 10~ 3 , 
much weaker than (f23|). 



2.4 Can QCD solve the strong CP problem itself ? 

Thus, theorists' task is to try to find a mechanism which would make CP con- 
servation in strong interactions natural. Two alternative approaches are logically 
possible. One can invoke physics beyond QCD (this approach will be discussed in 
the bulk of this review) or one can try a minimalistic standpoint and ask whether 
QCD could solve the strong CP problem itself, with no new physics. 

An obvious solution of the latter type exists: were one of the quarks massless, 
e.g., m u = 0, then all 9 effects would be unobservable. In this case there is a global 
Ua(1) symmetry of the chiral rotations of the u quark field, Ur — > exp(i/3) ur , ul — > 
exp(—ij3)uL, which can eliminate 9 altogether. However, m u = does not go 



through on phenomenological grounds and at present this scenario may be 

safely discarded. 

A more intricate solution could exist if confinement itself were to ensure the 
effects of the 9 term to be screened. As far as we know, this question was first 
raised by A. Polyakov shortly after the discovery of the strong CP problem. His 
argumentation was as follows. The 9 term in the action is the integral over the 
full derivative, which can be operative only if there are long-range components of 
the gauge fields. In the quasiclassical approximation such components are certainly 
present, as is evident from instanton calculations. However, this approximation 
misses the most salient feature of QCD, color confinement, which might eliminate 
long-range interactions (i.e. "screen" color) and make the integral over the full 
derivative vanish. That's exactly what happens in the (l+2)-dimensional Polyakov 
model of color confinement [19| 



This issue was studied in Ref. H17| , with the negative conclusion. The problem 
is that there are two effects in QCD which depend on the above "screening": the 
heaviness of the rf mass (the so-called U(l) problem) and CP conservation/violation. 
Even though we do not know precisely how exactly the confinement mechanism 
works in QCD, we know for a fact that rf is split from the octet of the Goldstone 
bosons. This knowledge (plus some reasonable arguments regarding the value of 
chiral and/or 1/N C corrections) is sufficient to show that confinement in QCD does 
not eliminate the 9 dependence and thus does not solve the strong CP problem. We 
briefly outline this line of reasoning below. 

The flavor singlet meson, the rj', is significantly heavier than the flavor octet 
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Goldtones, m n i ~ 958 MeV. As was shown by Weinberg in the pre-QCD era, were 
the rj' a Goldstone, its mass would be constrained by m v i < ^/Sm^. This suggests 
that, unlike the octet of the Goldstone bosons, the rj is not massless in the chiral 
limit (unless chiral expansion is invalid). This is the only fact we will need. 

An extra contribution to the rj' mass comes from nonperturbative effects due to 
the axial anomaly, as was exemplified |2(J by the same instantons. 

To quantify the effect on the theoretical side, let us introduce the correlator of 
the topological charge densities (in the literature it is referred to as the topological 
susceptibility) 

X = -i J d 4 x{0\TQ(x)Q(0) |0) , (24) 



where 



Q = 32 ^2 G% G^u ■ (25) 



According to the low-energy theorem derived in Refs. |2T|, [22|] in the leading 
approximation of the 1/N C expansion, the quantity X is saturated by the rf contri- 
bution implying the following formula for the rj' mass: 

K = ^ + ^ + {-^) > ^ 

where the topological susceptibility X on the right-hand side is evaluated in pure 
gluodynamics, the Yang-Mills theory with no light quarks. In order for the rj' mass 
to be nonzero in the chiral limit, X should be nonzero in pure gluodynamics. 

A substantial amount of theoretical evidence is accumulated in the last 20 years 
showing that the topological susceptibility in pure Yang-Mills does not vanish. The 
lattice m and the QCD sum rule studies g| yield X ~ (180 MeV) 4 ^ 0. This 



successfully takes care of the U(l) problem. 

Having nonzero topological susceptibility in pure gluodynamics means, per se, 
that there is a sensitivity to the parameter 9 in this theory. Indeed, X is nothing 
but the second derivative of the vacuum energy with respect to 9 taken at 9 = 



d 2 flnZo 

A, — 



d9 2 \ V 



(27) 

9=0 



In the theory with the light quarks included, the topological susceptibility can be 
calculated by applying the chiral perturbation theory (see e.g. Hi~7|). As expected 



on general grounds, in this case X oc m q provided the rj' is split from the octet 
of the Goldstones; otherwise the rf contribution cancels the 0(m q ) term in the 
topological susceptibility. However, in the theory with the light quarks it is much 
more instructive to calculate directly CP odd decay rates, for instance the rate of 
r] — ► 7r + 7r~. This amplitude is forbidden by CP. In the same way as with the 
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topological susceptibility, the chiral low-energy theorem yields |T7| a nonvanishing 
amplitude 0(m g ) provided the rj' is split from the octet of the Goldstones. 

Therefore, since QCD does solve the U(l) problem — it does split the rj' from 
the octet Goldstones — it cannot solve the strong CP problem |[l7j without help 
from outside f\. 



3 In search of a solution beyond QCD 
3.1 Peccei-Quinn mechanism 

The first dynamical mechanism solving the strong CP problem was proposed by 
Peccei and Quinn [571] . The main observation of Ref . p?| is as follows: if there is a 
U(l) axial symmetry in the theory 

q L ^e ia q L , qR^e^qjt, (28) 

then the 9 term can be removed from the Lagrangian, much like in the case of 
one massless quark discussed in the previous section. Below this symmetry will be 
referred to as U(1)pq. 

To see whether this symmetry is present in the Standard Model with one Higgs 
doublet let us consider the Yukawa sector and restrict ourselves to the first generation 
quarks (consideration in the general case is quite similar), 

XuQl^ur + \ d Q L <P c d R + H.c. + V(<f>+<f>), (29) 



2 The situation seems to be quite clear in this respect, nevertheless, an attempts to develop 
models of confinement that would solve the strong CP problem are not abandoned. Let us comment 
on a proposal of Ref. where two distinct "topological susceptibilities" are defined: the local 
and global ones. Let V c denote the volume at which the confinement effects take place, and V be 
the total volume of space time V c <C V — > oo. According to J25j, the local topological susceptibility 
is X loc = j Vc d A x (0\TQ(x) Q(0) |0) , while the global one is X gloh = J v d 4 x (0\TQ(x) Q(0) |0) . 
The author claims that the solution of the U(l) problem requires that X\ oc ^ 0, while the solution 
of the strong CP implies X g \ Q \> — 0, and both conditions may be dynamically compatible, so that 
both the strong CP and U(l) problems could be solved simultaneously. The underlying dynamics 
outlined in ]||| is a special interaction between instantons which "screens" them. 

There are a number of objections to this suggestions. First, it is the global topological sus- 
ceptibility that enters in the Witten-Veneziano relation and determines the rf mass modulo 1/N C 
corrections. Even if we forget about theoretical calculations of this quantity demonstrating that 
it does not vanish, we know that the rj' is split from the Goldstone octet because the Weinberg 
relation m r) < < y/Srn^ is grossly violated. (Of course, the validity of the chiral expansion is as- 
sumed; otherwise the Weinberg relation is meaningless.) If the rj' is split, there is no way out |jT?| : 
6>-induced effects are observable. 

In terms of the model suggested in J25| this means that if one were able to complete the calcu- 
lations at the level of physical observables, one would find that the rf is not split from the octet of 
the Goldstones or, more likely, that the required interaction between instantons is not sustainable. 
The latter variant was advocated in 1291 . 
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where Ql is the left-handed SXJ(2) W quark doublet, 



is the charge 



conjugate Higgs field, and V denotes the Higgs potential. Although the first term 
in this expression is invariant under the transformations 



Ql 



e ia q L 



Qr 



e~ ia q R 



(30) 



there is a second term which is not invariant under (RO) since <p c transforms as 
conjugate to <fi. Therefore, the one Higgs doublet SM is not invariant under U(1)pq 
and the strong CP cannot be solved. 



However, as was pointed out in Ref. IpTfl , the required U(1)pq symmetry is 
present in SM with two Higgs doublets — let us call them <fi and x- in this case the 
Yukawa sector for the first generation quarks reads 

KQl<I>u r + Xd Ql X* dn + H.c. + X + X, (0 + X) (x + <t>)) ■ (31) 

It is invariant under 



Ql 



e ia q L 



Qr 



e Qr 



X 



e~ 2ia X- 



(32) 



This fact can be used to solve the strong CP problem [2~7|. The symmetry ( |3"2"D is 
explicitly broken by the axial anomaly. As a result, instantons induce an effective 
potential for the 9 term. The potential can be calculated in a certain approximation 
] . The crucial model- independent fact is that the resulting potential is minimized 



at a zero value of the CP violating phase [27 . 



3.2 Weinberg-Wilczek axion 

When the Higgs fields develop vacuum expectation values (VEV's) the local elec- 
troweak symmetry group is spontaneously broken. This gives masses to the interme- 
diate W ± and Z vector bosons. Simultaneously, the global [7(1)pq is spontaneously 
broken too. Spontaneous breaking of the global symmetry leads to the emergence 
of a massless Goldstone boson, the axion in the present case [0, 0]. In SM with two 
Higgs doublets the axion is given by the following superposition 

a = - (Vm0o - v x lmxo ) , (33) 

v 

where 0o and xo denote the neutral components of the Higgs doublets. Moreover, 
v = yjvj, + v^, ~ 250 GeV, and and v x are the vacuum expectation values of 
and X-, respectively. In this approximation the axion is massless. However, as we 
mentioned above, nonperturbative QCD effects (such as instantons) give rise to a 
potential for the axion. Hence, the axion acquires a nonzero mass which can be 
estimated as follows [l], ^ 

T Til 

m a ~ ~ 100 KeV. (34) 

v 
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Moreover, the axion decay constant is 1/v. Therefore, we see that the Weinberg- 
Wilczek (WW) axion mass and decay constant are tied to the electroweak symmetry 
breaking scale v. This turns out to be too much of a constraint, and, as we will 
discuss in Sec. 3.5, the WW axion is excluded on the basis of existing experimental 
data. 



3.3 KSVZ axion 

If the scale of PQ symmetry breaking is much higher than the electroweak scale, 
then according to (|34|), the axion is much lighter and its decay constant is much 



smaller. Such an "invisible" axion would not be in conflict with experimental data. 
A scenario with the harmless axion was first proposed in Refs. [^SJ and |I7| (the 



KSVZ axion) . In the latter paper it was called phantom axion. Needless to say that 
to untie the axion from the electroweak scale one has to decouple the corresponding 
scalar fields from the known quarks and couple them to hypothetical (very) heavy 
fermion fields carrying color. 

In more detail, one introduces a complex scalar field $ coupled to a hypothetical 
electroweak singlet, a quark field Q in the fundamental representation of color SU(3), 

AC = <$>Q R Q L + H.c. . (35) 

The modulus of <£> is assumed to develop a large vacuum expectation value f/y/2, 
while the argument of $ becomes the axion field a, modulo normalization, 

a(x) = fa(x) , a(x) = Arg$(x) , / » A . (36) 

Then the low-energy coupling of the axion to the gluon field is 



AC = ja^- 2 G%G% , (37) 



so that the QCD Lagrangian depends on the combination 9 + a(x). 

In general, one could introduce more than one fundamental field Q, or introduce 
them in a higher representation of the color group. Then, the axion-gluon coupling 



(p7|) acquires an integer multiplier N, 



AC' = jaN^- 2 G%G%. (3f 



This factor N (not to be confused with the number of colors iV c nor with M of 
extended super symmetry) is sometimes referred to as the axion index. The minimal 
axion corresponds to N — 1. In the general case the QCD Lagrangian depends on 
the combination 9 + Na(x). As previously, nonperturbative QCD effects generate 
a potential for 9 + Na(x). The latter is minimizes at the value 9 + iVo; vac = 0, i.e., 
the strong CP problem is automatically solved. 
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3.4 ZDFS axion 



An alternative way to introduce an "invisible" axion was proposed in Refs. p9l and 



, (the ZDFS axion). In this proposal one maintains the PQ symmetry of the two 



doublet SM but separates the scales of the PQ and electroweak breaking G9, 30 



To this end the SM Lagrangian is extended - a scalar SM singlet field S is added, 
\ u Q L (pu R + X d QLX*d R + H.c. + 

V(<P + , X + X , (0 + X) , S + S , £ 2 ) . (39) 

One notes that this expression is invariant under the following axial transformations 

q L ^e ia q L , q R ^ e - ia q R , 0^e 2ia 0, x ^e" 2lQ X , £ - e 2lQ S . (40) 

Upon spontaneous breaking of this symmetry the Goldstone particle, an axion, 
emerges as a superposition 

a = — (f^Im^o - v x lmxo + ImS) , (41) 



where V = Jv^ + + t>|, and v^, v x and t>s are the vacuum expectation values 
of <p, x an d S, respectively. The vacuum expectation value of £ does not have to 
be related to the electroweak symmetry breaking scale. In fact, it can be as large 
as the GUT scale. If so, the axion is light and its decay constant is tiny. We will 
discuss experimental bounds on these quantities in the next section. 

3.5 Constraints on the axion mass 

As we discussed in the previous sections, the PQ symmetry is explicitly broken by the 
axial anomaly. Therefore, the axion is a pseudo Goldstone boson. Nonperturbative 
QCD effects induce the axion mass. For further discussions it is convenient to 
parametrize the axion mass as follows 

m a ~ 0.6eV 1Q7 ^ eV , (42) 

where / is the axion decay constant determined by the PQ breaking scale. 

In general, while discussing phenomenological constrains on the axion mass, one 
should distinguish between the KSVZ and ZDFS cases. The axion couplings to 
matter are different in these two scenarios. In particular, the KSVZ axion has no 
tree-level couplings to the standard model quarks and leptons. However, the aim 
of the present section is to summarize briefly an order of magnitude constraints on 
the axion mass. For this goal the effects which distinguish between the KSVZ and 
ZDFS axions will not be important (for detailed studies see Ref. [[H]] and citations 
therein). 
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The quantity 1/ / sets the strength of the axion coupling. Light axions can be 
produced in stars and a part of the energy of a star can be carried away by those 
axions. Stars can loose energy due to the production of light axions in the following 
possible processes 

(i) Nucleon-nucleon bremsstrahlung: N + N — > N + N + a; 

(ii) The Primakoff process: 7 <-» a conversion in the electromagnetic field of a 
nucleus; 

(iii) Photoproduction on an electron: 7 + e~ —>■ e~ + a ; 

(iv) Electron bremsstrahlung on a nucleus: e~ + (A, Z) — > e~ + a + (A, Z) ; 

(v) Photon fusion: 7 + 7^0; 

Unless 1/ f is really small, the emission of axions that are produced in the above 
reactions would lead to unacceptable energy loss by the star. This leads to the 
following lower bound on the axion decay constant / > 10 9 GeV. 

It is remarkable that cosmology puts an upper bound on / [021 K3L The 
latter comes about as follows. If / is too large then the axion coupling 1/ f is very 
small. As a result, during the course of cosmological evolution of the universe axions 
decouple early and begin to oscillate coherently. There are two major mechanisms 
by which the energy density stored in these oscillations can be dissipated - the 
Hubble expansion of the universe and the particle production by axions. However, 
if / > 10 12 GeV, neither of these mechanisms are effective (the axion coupling is 
too small). As a result, at some point of the evolution the axion energy density 
exceeds the critical energy density and over-closes the universe. In order for this not 
to happen one should impose the constraint / ^ 10 12 GeV. 

Summarizing, we obtain the following order of magnitude bounds on / and m a : 

10 9 GeV £ / £ 10 12 GeV, 1(T 6 eV £ m a £ 1(T 3 eV . (43) 
For further details see, e.g., Ref. |3"T|j . 



4 The vacuum structure in large N c gluodynamics 



The early studies |35]| of the chiral Ward identities in QCD revealed that the vacuum 



energy density depends on the vacuum angle 9 through the ratio 9/Nf, where Nf 
is the number of quarks with mass m q <C A. Shortly after it was shown in Refs. 
]36] and [|37| that this structure occurs naturally, provided that there exist Nf states 
in the theory such that one of them is the true vacuum, while others are local 
extrema; all are intertwined in the process of "the 9 evolution." Namely, in passage 
from 9 = to 9 = 2tt, from 9 = 2ir to 9 = Air, and so on, the roles of the above 
states interchange: one of the local extrema becomes the global minimum and vice 
versa. This would imply, with necessity, that at 9 = kir (where k is an odd integer) 
there are two degenerate vacuum states. Such a group of intertwined states will be 
referred to as the "vacuum family." The crossover at 9 = tt, 37r, etc. is called the 



D ashen phenomenon p8 |. 
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This picture was confirmed by a detailed examination of effective chiral La- 
grangians |36], |37|, |39], [10| (for a recent update see [§4]]). For two and three light 



quarks with equal masses it was found that the vacuum family consists of two or 
three states respectively; one of them is a global minimum of the potential, while 
others are local extrema.Q At 9 = n the levels intersect. Thus, Crewther's depen- 



dence |35| on 8/Nf emerges. 

On the other hand, the examination of the effective chiral Lagrangian with the 
realistic values of the quark masses, m d /m u ~ 1.8 , m s /m d ~ 20, yields |$6], [J7|, |3J 
a drastically different picture - the vacuum family disappears (shrinks to one state); 
the crossover phenomenon at 9 = n is gone as well. 

This issue remained in a dormant state for some time. Recently arguments were 
given that the "quasivacua" (i.e. local minima of the energy functional), which to- 
gether with the true vacuum form a vacuum family, is an indispensable feature of 



gluodynamics. The first argument in favor of this picture derives |E| from super- 
symmetric gluodynamics, with supersymmetry softly broken by a gluino mass term. 
The same conclusion was reached in Ref. f43[ based on a D-brane construction in 
the limit of large N c . One can see that in both approaches the number of states in 
the vacuum family scales as N c . In fact, in Ref. the expression for the theta 
dependence of the vacuum energy in the large N c pure Yang-Mills (YM) model was 
derived from a D-brane construction. It has the following form 



8(6) = C min fc (0 + 2nk) 2 + O (J^j , (44) 

where C is some constant independent of N c and k stands for an integer number. 
This expression has a number of interesting features which might seem a bit puzzling 
from the field theory point of view. Indeed, in the large N c limit there are N% degrees 
of freedom in gluodynamics, thus, naively, one would expect that the vacuum energy 
density in this theory scales as ~ N%. However, the leading term in Eq. (|44|) scales 
as ~ 1. As a natural explanation, one could conjecture that there should be a 
colorless massless excitation which saturates the expression for the vacuum energy 
density fl44]) . However, pure gluodynamics generates a mass gap and there are no 
physical massless excitations in the model. Thus, the origin of Eq. ([14]) seems to be 
a conundrum. We will discuss and resolve this puzzle in the next section. 

Note, that an additional argument in favor of the vacuum family may be found 



in a cusp structure which develops once one sums up j44j sub-leading in 1/N C terms 
in the effective r/ Lagrangian. At large N c = oo the states from the vacuum family 
are stable, and so are the domain walls interpolating between them j|3], [15]. 

When N c < oo the degeneracy and the vacuum stability is gone, strictly speak- 
ing. It is natural to ask what happens if one switches on the axion field. This 



3 We stress that the states from the vacuum family need not necessarily lie at the minima of 
the energy functional. As was shown by Smilga [[ll| , at certain values of 9 some may be maxima. 
Those which intersect at 9 = kit (k odd) are certainly the minima at least in the vicinity of 9 = kit. 

4 It has been conjectured long time ago in |3q|. 
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generically leads to the formation of the axion domain walls. The axion domain 
wall HU presents an excellent set-up for studying the properties of the QCD vac- 
uum under the 9 evolution. Indeed, inside the axion wall, the axion field (which, in 
fact, coincides with an effective 6) changes slowly from zero to 2n. The character- 
istic length scale, determined by the inverse axion mass m." 1 , is huge in the units 
of QCD, A -1 . Therefore, by visualizing a set of spatial slices parallel to the axion 
wall, separated by distances ^> A -1 , one obtains a chain of QCD laboratories with 
distinct values of 9 e & slowly varying from one slice to another. In the middle of the 
wall 8 cS = 7r. 

Intuitively, it seems clear that in the middle of the axion wall, the effective value 
of 9 e ff = 7r. Thus, in the central part of the wall the hadronic sector is effectively in 
the regime with two degenerate vacua, which entails a stable gluonic wall as a core 
of the axion wall. In fact, we deal here with an axion wall "sandwich." Its core is 
the so-called D-wall, see P7fl . 

Below we will discuss this idea more thoroughly. We also address the question 
whether this phenomenon persists in the theory with light quarks, i.e., in real QCD. 
Certainly, in the limit N c = oo the presence of quarks is unimportant, and the 
axion wall will continue to contain the D-wall core. As we lower the number of 
colors, however, below some critical number it is inevitable that the regime must 
change, the gluonic core must disappear as a result of the absence of the crossover. 
The parameter governing the change of the regimes is A/N c as compared to the 
quark mass m q . At m q <C A/N c , even if one forces the axion field to form a wall, 
effectively it is screened by a dynamical phase whose origin can be traced to the rj', 
so that in the central part of the axion wall the hadronic sector does not develop 
two degenerate vacua. The D-walls cannot be accessed in this case via the axion 
wall. 

The issue of hadronic components of the axion wall in the context of a potential 
with cusps |Q] were discussed in p8| , |4"9| , |50|| . However, the gluonic component of 



the axion walls was not studied. The rf component in the axion walls was discussed 
in Ell 031. 



4.1 Arguments from supersymmetric gluodynamics 

First we will summarize arguments in favor of the existence of a nontrivial vacuum 
family in pure gluodynamics. 

The first indication that the crossover phenomenon may exist in gluodynamics 
comes |42j from supersymmetric Yang-Mills theory, with supersymmetry being bro- 



ken by a gluino mass term. The same conclusion was reached in Ref. |43[ based on 
a D-brane construction in the limit of large N c . In both approaches the number of 
states in the vacuum family is N c . 

The Lagrangian of softly broken supersymmetric gluodynamics is 



14 



+ 22^ G %G%>> ( 45 ) 

where m is the gluino mass which is assumed to be small, m< A (here we rescaled 
the gluon and gluino fields so that 1/g 2 appears as a common multiplier in the 
Lagrangian) . 

There are N c distinct chirally asymmetric vacua, which (in the m = limit) are 
labeled by 

(A 2 ), = iV c A 3 expfz^^V £ = 0,1,...,N C -1. (46) 



N, 



c 



At m = there are stable domain walls interpolating between them |52||. Setting 



m 7^ we eliminate the vacuum degeneracy. To first order in m the vacuum energy 
density in this theory is 

E™ (A 2 ) + H.c. = -m iV c 2 A 3 cos . (47) 

g 2 N c 

Degeneracy of the vacua is gone. As a result, all the metastable vacua will decay 
very quickly. Domain walls between them, will be moving toward infinity because 
of the finite energy gradient between two adjacent vacua. Eventually one ends up 
with a single true vacuum state in the whole space. 

For each given value of 9 the ground state energy is given by 

8(9) = min £ j —m N 2 A 3 cos ^^^ j . (48) 

At 9 = 7r, 37r, we observe the vacuum degeneracy and the crossover phenomenon. 
If there is no phase transition in m, this structure will survive, qualitatively, even 
at large m when the gluinos disappear from the spectrum, and we recover pure 
gluodynamics. 



Based on a D-brane construction Witten showed [|3] that in pure SU(Af c ) (non- 
supersymmetric) gluodynamics in the limit A^ c — » oo a vacuum family does exist :0 
the theory has an infinite group of states (one is the true vacuum, others are non- 
degenerate metastable "vacua") which are intertwined as 9 changes by 2irx (integer), 
with a crossover at 9 = ttx (odd integer). The energy density of the k-th state from 
the family is 

E k {9) = N 2 C A 4 F ( 9 -^) , (49) 



5 This was shown in Ref. [(43| assuming that there is no phase transition in a certain parameter 
of the corresponding D-brane construction. In terms of gauge theory, this assumption amounts of 
saying that there is no phase transition as one interpolates to the strong coupling constant regime. 
Thus, the arguments of |0 have the same disadvantage as those of SUSY gluodynamics where 
one had to assume the absence of the phase transition in the gluino mass. 
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where F is some 27r-periodic function, and the truly stable vacuum for each 9 is 
obtained by minimizing 8k with respect to k, 



8(9) =AT c 2 A 4 min fc F^±^) , (50) 



v c 



much in the same way as in Eq. 

At very large N c Eq. (|50|) takes the form 



8(9) = A 4 min fc (9 + 2nkj 2 +Oi — \ . (51) 



N, 



The energy density 8(9) has its absolute minimum at 9 = 0. At N c — oo the 
"vacua" belonging to the vacuum family are stable but non-degenerate. To see that 
the lifetime of the metastable "vacuum" goes to infinity in the large N c limit one can 
consider the domain walls which separate these vacua p5|, |53f| . These walls are seen 



as wrapped D-branes in the construction of [43], and they indeed resemble many 
properties of the QCD D-branes on which a QCD string could end. We refer to them 
as D-walls because of their striking similarity to D2-branes. The consideration of 
D-walls has been carried out |45| and leads to the conclusion that the lifetimes of 



the quasivacua go to infinity as exp (const N^). 

Moreover, it was argued [^, f|7| that the width of these wall scales as 1/N C both, 



in SUSY and pure gluodynamics. To reconcile this observation with the fact that 



masses of the glueball mesons scale as N®, we argued |47| that there should exist 
heavy (glue) states with masses oc N c out of which the walls are built. The D-brane 
analysis |55| , effective Lagrangian arguments and analysis of the wall junctions [56j, 



support this interpretation. These heavy states resemble properties of DO-branes. 
The analogy is striking, as DO-branes make D2-branes from the standpoint of the 
M(atrix) theory [p7B , so these QCD "zero-branes" make QCD D2-branes (domain 
walls) .Q The distinct vacua from the vacuum family differ from each other by a 
restructuring of these heavy degrees of freedom. They are essentially decoupled 
from the glueballs in the large N c limit. 
Now we switch on the axion 

AC = l -f(d»a)(d»a) + ^ G%G% , (52) 



with the purpose of studying the axion walls. The potential energy 8(9) in Eq. ( |50|) 
or ( |5"TD is replaced by 8(9 + a). 

Since the hadronic sector exhibits a nontrivial vacuum family and the crossover^ 
at 9 = 7r, 37r, etc., strictly speaking, it is impossible to integrate out completely the 
hadronic degrees of freedom in studying the axion walls. If we want to resolve the 
cusp, near the cusp we have to deal with the axion field plus those hadronic degrees 



6 See also closely related discussions in Ref. |p8| . 

7 For nonminimal axions, with N > 2, the crossover occurs at a = kw/N. 
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of freedom which restructure. In the middle of the wall, at a = ir, it is mandatory 
to jump from one hadronic vacuum to another - only then the energy of the overall 
field configuration will be minimized and the wall be stable. Thus, in gluodynamics 
the axion wall acquires a D-wall core by necessity. 

One can still integrate out the heavy degrees of freedom everywhere except a 
narrow strip (of a hadronic size) near the middle of the wall. Assume for simplicity 
that there are two states in the hadronic family. Then the low-energy effective 
Lagrangian for the axion field takes the form (|83|). The domain wall profile will also 
exhibit a cusp in the second derivative. The wall solution takes the form: 

{8 arctan (e ma2 tanf ) , at z < 
t , v (53) 

-2tt + 8 arctan (e maZ tan^J , at z > , 

at N c = 2 (the wall center is at z — 0). 

Examining this cusp with an appropriately high resolution one would observe 
that it is smoothed on the hadronic scale, where the hadronic component of the 
axion wall "sandwich" would become visible. The cusp carries a finite contribution 
to the wall tension which cannot be calculated in the low-energy approximation but 
can be readily estimated, T core ~ A 3 iV c . In subsequent sections we will examine this 
core manifestly in a toy solvable model. 



4.2 Argumens from D-brane construction 



The theta dependent vacuum energy (53) is related to the correlator measuring the 
vacuum fluctuations of the topological charge. The question which arises here is 
whether this can be seen from the original string theory computation [53|. We are 
going to discuss below how the string theory calculation suggests that the vacuum 
energy (03J) should indeed be related to the vacuum fluctuations of the topological 
charge. In fact, we argue that this is related to the instantons carrying .DO-brane 
charge in the Type IIA fourbrane construction of the four- dimensional YM model. 

In general, a great deal of information can be learned on nonperturbative phe- 
nomena in four- dimensional gauge theories by obtaining these models as a low-energy 
realization of certain D-brane configurations flSSH , and/or using a duality of large N 
superconformal gauge theories and string theory compactified on certain spaces (see 
Refs. |6(| and |61], ^]). This duality, being a powerful technique, has also been gen- 
eralized for the case of non-supersymmetric models p3| |. This was applied to study 
various dynamical issues in large N c pure Yang- Mills theory |53], |55L EBL |57|, |5SL . 

To begin with let us recall how the theta dependent vacuum energy appears in 
the brane construction of the four-dimensional YM model [E3I. One starts with 



Type IIA superstring theory on Ai = R 4 xS'x R 5 , with N c coincident D4-branes 



|63| . The ZM-brane worldvolume is assumed to be R 4 x S 1 and the fermion bound- 
ary conditions on S 1 are chosen in such a way that the low-energy theory on the 
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worldvolume is pure non-supersymmetric U(N C ) YM theory ||63|| . In the dual de- 
scription, the large N c limit of the SU(N C ) part of this theory can be studied by 
string theory on a certain background |BD|, ^T|, ^3]. It was shown in Ref. |I3j] that 
the theta dependent vacuum energy ([0]) arises in the dual string description due to 
the U(l) gauge field Bm, M = 1, .., 5. To find out what this corresponds to in the 
original gauge theory language recall that this U(l) field is nothing but the Ramond- 
Ramond (RR) one-form of Type IIA theory. Furthermore, once the gauge theory is 
realized in the Type IIA fourbrane construction, the Wess-Zumino-Witten (WZW) 
term present in the worldvolume effective action defines the correspondence between 
the gauge theory operators on one side and the string theory Ramond-Ramond fields 
on the other side. In the case at hand the worldvolume WZW term looks as follows: 



'WZW : 



tt 2 Jn 



BATrGAG, (54) 



where Q denotes the worldvolume of a wrapped fourbrane, Q = R 4 x S 1 . In 
accordance with the general principles of the large N c AdS/CFT correspondence 



S2| the classical action for the RR one-form on the string theory side defines 



the YM correlation functions of the composite operator GG (since this is the opera- 



tor which couples to the corresponding RR field in fl54)). Thus, it is not surprising 
that the theta dependent vacuum energy which is defined by the RR one-form in 
the string theory calculation is related to the nonzero value of the topological sus- 
ceptibility in the gauge theory studies. The physical reason for this correspondence, 
as we have mentioned above, is the special property of the gauge theory instantons 
in the fourbrane construction. Indeed, in accordance with ( CT ) the RR one-form 
couples to the topological charge density GG, on the other hand the RR one-form 
couples by definition to DO-branes. Thus, the gauge theory instantons in this case 
carry zerobrane charge. This is the physical reason for the correspondence discussed 
above. 



4.3 Derivation of vacuum energy in QCD 

The aim of this section is to derive Eq. ((44|) in pure YM model and, in particular, 
to identify the degrees of freedom which are responsible for the theta dependent 
vacuum energy density. 

In the quasi-classical approach the theta dependence can be calculated using 
instantons 0. In a simplest approximation of non-interacting instantons the theta 
angle enters the Euclidean space partition function in the following form: 

exp(-^±i0) =exp(-AT c ^±^) , (55) 

where g stands for the strong coupling constant. A denotes the 't Hooft's coupling 
X = N c g 2 which is kept fixed in the large N c limit. The expression above vanishes 
in the large N c limit, so does the theta dependence in (|55"D . However, this conclusion 
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cannot be extrapolated to the infrared region of the model. The limitations of the 
expression fl5"5|) prevent one to do so. Indeed, the quasi-classical approximation is 
valid in the limit of small coupling constant (see, for instance, discussions in Ref. 
|ffO|). Once quantum corrections are taken into account the coupling constant g 2 in 
fl55|) becomes a scale dependent quantity. In fact, it will depend upon an instanton 
size g 2 = g 2 {p)- For small size instantons the running coupling is small and the 
quasi-classical approximation in ( |55|) holds. However, for large size instantons, i.e. 
large couplings, it is not even clear whether the notion of a single instanton is a 
legitimate approximation. The overlap between instantons can be big in this case 
and some more complicated field configurations should be relevant for the description 
of physical phenomena [fflH . In any event, the expression fl55|) is no longer reliable in 
the strong coupling limit. Thus, the conclusion that the theta dependence goes away 
in the large N c limit cannot be justified. One way to study the infrared region is 
to look for some appropriate composite colorless excitations for which the notion of 
an asymptotic state can be used. We will start by searching for these excitations in 
pure Yang-Mills theory. To proceed, let us recall that the topological susceptibility, 
X, is a nonzero number in pure gluodynamics (we rewrite it in the following form) 

X = -i yV 9"<0|T K„(x)K v (0)\0) d 4 x ^ . (56) 

Here, as before denotes the Chern-Simons current. As we discussed earlier, the 
value of X in large N c pure YM theory determines the rf meson mass in full QCD 
with massless quarks via the Witten-Veneziano formula m 2 ,f 2 , oc X, with f v i being 
the 7/ meson decay constant |21, p2|f|. 



In what follows it will prove convenient to introduce a new variable by rewriting 
the expression for the topological charge density Q in terms of a four-index (four- 
form) tensor field H^ 1 " 1 ^: 

Q = E ^ Va " , (57) 
where the four-form field H^ ua ^ is the field strength for the three-form potential 

n 

= d^C ua/ 3 — d u C^ a/ 3 — d a C ufl/ 3 — dpCva^. (58) 

The C^a field is defined as a composite operator of the gluon fields A^: 

C, ua = ^- 2 {Ap v Al - A a v d,Al - A a a d v Al + 2f abc A»A b u Al). (59) 

Here, f a t, c denote the structure constants of the corresponding SU (N c ) gauge group. 
The right-left derivative in this expression is defined as AdB = A(dB) — (dA)B. 

8 Below, unless otherwise stated, we will not distinguish between X and its large N c limit. The 
constant contact term in the definition of X will also be omitted for simplicity. 
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Notice, that the C va p field is not a gauge invariant quantity; if the gauge transfor- 
mation parameter is A a , the three-form field transforms as C va p — > C ua p + d v K a p — 
— dpA au , where A Q/3 oc A^dpA 0, — A a pd a k a . However, one can check that the 
expression for the field strength H^ ua p is gauge invariant. 

It has been known for some time J72J that the C va $ field propagates long-range 
correlations if the topological susceptibility is nonzero in the theory. The easiest way 
to see this is to turn to the notion of the Kogut-Susskind pole f73fl . Let us consider 
the correlator of the vacuum topological susceptibility at a nonzero momentum. In 
this case X is defined as the zero momentum limit of the correlator of two Chern- 
Simons currents multiplied by two momenta: 

X — —i lim <f<f / e iqx (0\TKJx)K u (0)\0)d 4 x. (60) 
g— >-0 J 

Since this expression is nonzero, it must be that the correlator of two Chern-Simons 
currents develops a pole as the momentum vanishes, the Kogut-Susskind pole ]73 



Given that the correlator of two Chern-Simons currents has a pole and that the 
Chern-Simons current and the three-form C va p field are related, one concludes that 
the C ua/ 3 field also has a nonzero Coulomb propagator |72| . Thus, the C va p field 



behaves as a massless collective excitation propagating a long-range interaction [72 
These properties, in the large N c limit, can be summarized in the following effective 
action for the C va p field: 

1 f 9 f 

Seff = ~ 2 4 , x J H* uaf3 d A x - - J ^ C vaP dx u A dx a A dx p + High dim. (61) 

The first term in this expression yields the correct Coulomb propagator for the 
three-form C va @ field. The second term is just the usual CP odd 9 term of the initial 
YM action written as a surface integral at spatial infinity dT. Notice that higher 
dimensional terms are not explicitly written in this expression. There might be two 
types of higher dimensional contributions in (|6T]) . First of all, there are terms with 
higher powers of derivatives of the fields. These terms are suppressed by momenta 
of the "massless" three-form field and do not contribute to the zero momentum 
vacuum energy of the system. In addition, there might be higher dimensional terms 
with no additional derivatives. In the next section we will present some of these 
contributions and show that they are suppressed by higher powers of 1/N C . 

In what follows we are going to study the large N c effective action given in Eq. 
(|Bl"l) |. In particular, we will calculate the ground state energy of the system in the 
large N c limit using the effective action (0). In fact, we will derive Eq. (|44|). 



9 The action (|6l]) is not an effective action in the Wilsonian sense. It is rather related to the 
generating functional of one-particle-irreducible diagrams of the composite field in the large A^ c 
limit. The effective action in Eq. (^lj) is not to be quantized and loop diagrams of that action are 
not to be taken into account in calculating higher order Green's functions. The analogous effective 



action for the CP even part of the theory was constructed in Refs. J74], 75 
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Before we turn to this calculation let us mention that Maxwell's equations for 
a free four-form field-strength H^ ua/3 yield only a constant solution in (3 + 1)- 
dimensional space-time [[7E| ]. The reason is as follows. A four-form potential has 
only one independent degree of freedom in four-dimensional space-time, let us call 
it E. Then, the four Maxwell's equations written in terms of the £ field require that 
this field is independent of the all four space-time coordinates, thus the solution 
can only be a space-time constant. As a result, the free H^ ua/3 field propagates no 
dynamical degrees of freedom in (3 + l)-dimensions. However, this field can be re- 
sponsible for a positive vacuum energy density in various models of Quantum Field 
Theory (see Ref. [[HJ). Thus, studying classical equations of motion for the H^ UQ/3 
field one can determine the value of the ground state energy given by configurations 
of H^ ua/3 . We are going to solve explicitly the classical equations of motion for the 
effective action (|61|) . Then, the energy density associated with those solutions will 
be calculated. 

Let us start with the equations of motion. Taking the variation of the action 
(|B"l~l) with respect to the C va p(z) field one gets 

VH^apiz) = e X J S (4 \z - x) dx u A dx a A dx p . (62) 



This equation can be solved exactly in four-dimensional space-time |[T6|| . The solu- 
tion is the sum of a particular solution of the inhomogeneous equation and a general 
solution of the corresponding homogeneous equation: 

Hfj, va p(z) = 6 X J 5 ( - 4 \z — x) dx^ A dx v A dx a A dxp + b e^ap. (63) 

The integration constant b, if nonzero, induces an additional CP violation beyond 
the existed 6 angle. However, periodicity of the 9 angle with respect to shifts by 
2tt x (integer) allows for some nonzero b proportional to 27rZ. As a result, the general 
solution to the equation of motion reads as follows: 

= -(0 + 2irk) X e^ af3 . (64) 

Thus, the different vacua are labeled by the integer k and the order parameter for 
these vacua in the large N c limit can be written as: 

(GG} k = (6 + 2nk) X. (65) 

As a next step let us compute the vacuum energy associated with the solution given 



in Eq. (0). The density of the energy-momentum tensor for the action (^1|) takes 
the form 



QfiV ~ 3! X 



H fia/3T H u aPT - -g^H^p^j . (66) 
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Using the expression flS^D one calculates the corresponding energy density^] 



£k(0) = ^(6 + 2nk) 2 X. (67) 

Since the H^ ua p field does not propagate dynamical degrees of freedom the expression 
above is the total energy density of the system given by the action (|BTD 0. 

Before we go further let us stop here to discuss some of the consequences of Eq. 
(|67|). First of all, let us notice that the result fl67D , as well as Eq. (f44[) , is only valid in 
the limit of infinite N c . Below we will calculate subleading order corrections to Eqs. 



, |67j ) and argue that these expressions can also be used as a good approximation 



for large but finite N c . The constant C emerging in (M ) is related to the topological 
susceptibility as follows: 



Thus, the vacuum energy (|44l , |6~7D is defined by vacuum fluctuations of the topological 
charge measured by X . 

The crucial feature of (^7]) is that it defines an infinite number of vacua. The 
true vacuum is obtained by minimizing ( |6"TD with respect to k as in 

So{6) = -X min fe (0 + 2nk) 2 . 
This expression is periodic with respect to shifts of 9 by 27rZ and is also a smooth 



function of 9 except for 9 = n j|3| (see also discussions below). Thus, there are an 
infinite number of the false vacua in the theory [13 . The fate of these states will be 
discussed in section 5. 

Let us now consider full QCD with three quark flavors. We are going to write 
down a low-energy effective Lagrangian for this case and then gradually decouple 
quarks by taking the quark masses to infinity. The resulting effective Lagrangian 
should be giving the energy density for pure Yang-Mills theory. 

In the large N c expansion the effective Lagrangian of QCD with three flavors 
takes the form H79L B9l HUl : 



C(U, U*,Q) = C (U, U*) + ~i Q(x) Tr (in U - In U* 

' Q 2 (i) + 0Q(i) + 4= Tr (MU + M*U*) + . . . , (68) 



2 X ^ w ^ w 2V2 



10 Notice that the total YM energy density should contain some negative constant related to the 
nonzero value of the gluon condensate H]. This constant is subtracted from the expression for 
the energy discussed in this work. The energy density ( |67| ) is normalized as £0(8 = 0) = 0, and for 
k = was discussed in Q . 

11 One might wonder whether the same result is obtained if one treats 9 not as a constant 
multiplying Q in the Lagrangian, but as the phase that the states acquire under a topologically 
non-trivial gauge transformations. In this case the arbitrary integration constant in Eq. (|6^) has 
to be chosen in such a way which would guarantee a proper 8 dependence of the VEV of the 
topological charge density. This would leave the results of our discussion without modifications. 
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where U denotes the flavor group matrix of pseudoscalar mesons, Cq denotes the 
part of the Lagrangian which contains the meson fields only [fT9| , [39| , |40] , B is some 



constant related to the quark condensate, and M stands for the meson mass matrix 
(for recent discussions of the effective chiral Lagrangian approach see Ref. |80|). 
Higher order terms in (E5|) are suppressed by quadratic and higher powers of 1/N C . 
In order to study vacuum properties, we concentrate on the low-momentum approx- 
imation. The Lagrangian presented above can be used to solve the U(l) problem 



21, |22|]. Indeed, the field Q enters the Lagrangian in a quadratic approximation and 
can be integrated out. As a result, the flavor singlet meson, the r/, gets an additional 
contribution into its mass term. This leads to the Witten-Veneziano relation and 



the solution of the U(l) problem without instantons [21, 2^]. In the present case 
we would like to follow an opposite way. Namely, we are going to make quarks very 
heavy and integrate them out keeping the field Q in the Lagrangian. In the limit 
m q — ► oo one finds that M — > oo. Thus, the low-energy effective Lagrangian which 
is left after the mesons are integrated out will take the form: 

4 ff (Q) = ^Q 2 (x) + 0Q(x) + O(^, ±). (69) 

QCD c 

Rewriting the field Q in terms of the "massless" tensor C a ^ as in the previous 
section, one finds that the expression ( |69|) is nothing but the Lagrangian presented 
in (pip. Thus, the higher order terms neglected in fl6"T| ) which could contribute 
to the vacuum energy at zero momenta would correspond to higher corrections in 
l/N c . In fact, the subleading corrections to the effective Lagrangian fl6"B|) can also 
be found |8T|. These terms are proportional (with the corresponding dimensionful 
coefficients) to the following expressions: 

^ Q *TY (aW .), (70) 

The terms listed above are suppressed in the effective Lagrangian by the factor l/iV c 2 . 
As a next step, we can include the terms ( |7D| ) into the full effective Lagrangian and 
then integrate the heavy meson fields out. The net result of this procedure is that 
the terms proportional to Q 4 appear in the effective Lagrangian for pure YM theory. 
This, in its turn, modifies the equation of motion for the single component of H^ va p 
considered in the previous section. Performing the calculation of the vacuum energy 
in the same manner as discussed above we find the following result for the energy 
density: 

S k (6) = \x(0 + 2ntf + C -^*{e + 2nk) A + <D(±). (71) 

In this expression a constant emerges as a result of integration of the equation 
of motionF]. Notice that the topological susceptibility in the expression above is 



12 



The numerical value of this constant was recently calculated on the lattice |jS2| 
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also defined in the corresponding order in the large N c expansion: X = 2C + 
Xi/N c + X 2 /N^. Thus, the expressions ( |67l , |71~|) could in principle give a reasonable 
approximation for big enough but otherwise finite N c . The true vacuum energy 
density, £q(0), can be obtained by minimizing the expression (|7TD with respect to k 
as in (|4"4"[). Then, £q{6) satisfies the relation <9|£o(#) \e=o = X, no matter what is 
the value of the constant in (|7T|). 



5 Dynamics of false vacua 

In this section we will discuss the dynamics of the false vacua present in the theory. 
In accordance with ( p7| , [7T| ) there are an infinite number of vacua for any given value 
of the theta angle. Clearly, not all of these are degenerate. As we discussed, the true 
vacuum state is defined by minimizing the expressions (^O) with respect to k. All 
the other states are false vacua with greater values of the energy density. There 
is a potential barrier that separates a given false vacuum state from the true one. 
Thus, a false vacuum can in general decay into the true state through the process 
of bubble nucleation The decay rates for these vacua were evaluated in Ref. 

I5J. In this section we analyze the fate of the false vacua for different realizations 



of the initial conditions in which the system is placed. For the sake of simplicity 
we will be discussing transitions between the vacuum states labeled by k' and k for 
different values of these integers. The first two cases considered in this section were 
studied in Refs. [4~3] and the remaining of the section follows Ref. 



5.1 False vacua with k' ~ 1 

In this subsection we consider the system which in its initial state exists in a false 
vacuum with k' of order ~ 1. Let us start with the case when N c is a large but 
finite number so that the formula (^) (or (|7T|) ) is still a good approximation. Since 
there exists the true vacuum state with less energy, the false vacuum can "decay" 
into the true one via the bubble nucleation process. That is to say, there is a finite 
probability to form a bubble with the true vacuum state inside. The shell of the 
bubble is a domain wall which separates the false state from the true one. The 
dynamical question we discuss here is whether it is favorable energetically to create 
and expand such a bubble. Let us study the energy balance for the case at hand. 
While creating the shell of the bubble one looses the amount of energy equal to the 
surface area of the bubble multiplied by the tension of the shell. On the other hand, 
the true vacuum state is created inside the bubble, thus, one gains the amount of 
energy equal to the difference between the energies of the false and true states. The 



13 This decay can go through the Euclidean "bounce" solution |84 . Though the existence of 
the bounce for this case is not easy to understand within the field theory context, nevertheless, 
one could be motivated by the brane construction where this object appears as a sixbrane bubble 
wrapped on a certain space [ff3| . 
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energy balance between these two effects defines whether the bubble can be formed, 
and, whether the whole false vacuum can transform into the true one by expanding 
this bubble to infinity. Let us start with the volume energy density. The amount of 
the energy density which is gained by creating the bubble is[^]: 



AS = - X 

2 



(9 + 2vr) 



0' 



2vr X (6 + tt) 



Thus, AS scales as ~ 1 in the large N c limit as long as the volume of the bubble is 
finite. Let us now turn to the surface energy which is lost. This energy is defined 
as: 



E s = T D x (surface area) . (72) 

The tension of the wall between the adjacent vacua, Tjj, should scale as Tp ~ N c 
in the large N c limit. Hence, the surface energy will also scale as ~ N c . Thus, the 
process of creation of a finite volume bubble in the large N c limit is not energetically 
favorable. Indeed, the amount of energy which is lost while creating the shell is 
bigger than the amount which is gained. In terms of the false vacuum decay width 
this means that the width of this process is suppressed in the large N c limit |4"5" |: 



Volume 



oc exp 



(-«<), (73) 



where a stands for some positive constant E5[. Thus, one concludes that in the limit 



N r — > oo the false vacua with k' ~ 1 are stable [43, 45 . 



5.2 False vacua with k! ~ N c 

Here we study the fate of the false vacua with k! ~ iV c . We discuss a possibility 
that these vacua can decay into a state k with k! — k ~ N c and k' + k ~ A^ c . As 
in the previous subsection, we are going to study the energy balance for the bubble 
nucleation process. The amount of the volume energy density which is gained by 
creating such a bubble in the large N c limit scales as follows: 

AS oc X N 2 C . 

Thus, the volume energy which is gained increases as ~ N%. Let us now turn to 
the surface energy which is lost while nucleating a bubble. This is defined as E' s = 
T' D x (surface area), where T' D denotes the tension of the domain wall interpolating 
between the vacua labeled by k' and k. Since k' — k ~ N c these vacua are not 
neighboring ones. Thus, in general, there is no reason to expect that the tension 
of the wall interpolating between these vacua scales as ~ iV c . T' D might scale as 
~ N% at most (as the energy of a generic configuration in a model with N% degrees 

14 In this subsection we assume that ^ ±7r. The case 8 = tt will be considered below. 
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of freedom). However, even in the case when T' D ~ N% the volume energy which is 
gained is at least of the same order as the surface energy which is lost. Hence, it 
is energetically favorable to increase the radius of such a bubble (since the volume 
energy scales as a cubic power of the radius while the surface energy scales only as a 
quadratic power of the effective size). Thus, the bubble nucleation process will not 
be suppressed and the false vacua with k! ~ N c will eventually decay into the true 
ground state. Note, that the state k' = N c can as well decay into the neighboring 
vacuum k = N c — 1 which subsequently is allowed to turn into the ground state. 



5.3 Parallel domain walls 

In this subsection we consider the special case when all the vacua are present si- 
multaneously in the initial state of the model. This can be achieved, for instance, 
by placing in space an infinite number of parallel domain walls separating different 
vacua from each other. It is rather convenient to picture these walls as parallel 
planes. Each vacuum state is sandwiched between the corresponding two domain 
walls (two planes) separating this state from the neighboring vacua. Each domain 



is labeled by k and in accordance with ( |67| , |7TD is characterized by the corresponding 



value of the vacuum energy. Furthermore, the order parameter (GG) takes different 
values in these vacua in accordance with (|65|). Let us turn to the true vacuum state. 
For simplicity we assume that this state is given by k = (which corresponds to \9\ 
being less than 7r). The corresponding vacuum energy is the lowest one. Consider 
the two states which are adjacent to the true vacuum. These states have the energy 
density bigger than that of the true vacuum. Thus, there is a constant pressure 
acting on the domain walls separating the true vacuum from the adjacent false ones. 
This pressure will tend to expand the domain of the true vacuum. In fact, for large 
but finite N c , the pressure will indeed expand the spatial region of the true vac- 
uum by moving apart the centers of the domain walls sandwiching this state. The 
very same effect will be happening between any two adjacent vacua. Indeed, let us 
calculate the jump of the energy density between the two vacua labeled by k! and 
k: 

AS k , k = 2tt X (k' -k) (9 + n(k' + k)) . (74) 

As far as N c is large but finite, the walls will start to accelerate. Farther the wall 
is located from the true vacuum (i.e. larger the sum k' + k), bigger the initial 
acceleration of the wall is going to be; i.e., the walls will start to move apart from 
each other with the following initial acceleration: 

(k' -k) \6 + 7i(k' + k)] 
a k >k oc A YM M " ■ ( 75 ) 

For finite N c all the walls will be moving to spatial infinity and the whole space will 
eventually be filled with the true vacuum state. On the other hand, when N c — > oo 
the picture is a bit different. There are a number of interesting cases to consider: 
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First of all let us set k' — k = 1 and A;' , ~ 1. Then, in the limit iV c — > oo 
the acceleration — > 0. Thus, the neighboring walls stand still if they had no 
initial velocity. The physical reason of this behavior is as follows. Although there 
is a constant pressure of order ~ 1 acting on the wall, nevertheless, the wall cannot 
be moved because the mass per unit surface area of the wall tends to infinity in the 
limit N c — > oo. 

The second interesting case would be when the constant pressure produced by 
the energy jump between some neighboring vacua is of order ~ N c . In this case it 
will be possible to accelerate these walls up to the speed of light and send them to 
spatial infinity. Indeed, if k' — k = 1 but k' + k ~ N c , then the wall between these 
two vacua starts moving with a non-vanishing acceleration which scales as follows: 

ay k oc A YM 1 ^ - ~ 0{l). (76) 

Thus, these walls will eventually be approaching spatial infinity with a speed of light 
even in the limit of infinite N c . 

In addition to the effects emphasized above there might also be decays of the 
false vacua happening in each particular domain. As we discussed in the previous 
subsections, for large but finite N c all the false vacua will be nucleating bubbles with 
energetically favorable phases inside and expanding these bubbles to infinity. Thus, 
for large but finite N c , there are two effects which eliminate the false vacua: The 
moving walls are sweeping these states to infinity, and, in addition, these vacua are 
decaying via bubble nucleation processes. 

What happens for an infinite iV c ? As we learned above there are an infinite 
number of domains which will stay stable in that limit and the corresponding false 
vacua would not decay because of the exponential suppression. Thus, there are an 
infinite number of inequivalent spatial regions which are separated by domain walls. 
Consider one of the regions sandwiched between two domain walls. The three-form 
field C^a will couple to the walls and the large N c effective action for this case will 
look as follows: 

S = S cS + A*i / G nw dxl1 A dxV A dx", (77) 

i=k, fc+1 7w < 

where S e s is defined in (|EH]), /ij stands for the coupling of the three-form potential 
to a corresponding domain wall; Wj denotes the worldvolume of the wall. In this 
case the domain wall can be regarded as a source of the corresponding three-form 
potential. This is reminiscent to what happens in the large N c supersymmetric YM 
model |85| . 



5.4 Domain walls at 9 = it 

If 9 — 7r, the initial classical Lagrangian is CP invariant. Indeed, under CP trans- 
formations 9 = 7i goes into — ix. Since it and — ix angles are equivalent, CP is a 
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symmetry of the Lagrangian. However, in accordance with (BSI), this symmetry is 
spontaneously broken by the vacuum of the theory. Thus, one finds the following 
two degenerate true vacua: 

£k=o = £fc=-i = 2 % 7r 2 . (78) 

These two vacua are labeled by the order parameter (|65D. In the k = state 
(GG) = 7iX and in the k = — 1 state (GG) = —nX. As a result of the spontaneous 
breaking of a discrete symmetry there should be a domain wall separating these 
two vacua. Let us consider the case discussed in the previous section. Namely, 
let us choose the initial condition of the system as a state where all the possible 
vacua are simultaneously realized in space. That is, there are an infinite number 
of domain walls (parallel planes) dividing space into an infinite number of domains 
with different values of the vacuum energy density labeled by k. As we mentioned 
above, only two of these domains have equal minimal energy density given in (|7£|). 
The domain wall separating these two vacua, as we will see below, is somewhat 
special. In accordance with the discussions in the previous subsection, for large but 
finite N c , all the walls merging with the false vacua will tend to rush to spatial 
infinity. The final stable state of the model can be described as a space separated 
into two parts by a single domain wall. To the left (right) of the wall one discovers 
the phase with k — — 1 with the corresponding order parameter (GG) = —nX , 
and, to the right (left) of the wall one finds the state with k = and (GG) = %X. 
Vacuum energies of these two states are degenerate. 

In the case of infinite N c the picture is slightly different. As elucidated in the 
previous subsection, there will be an infinite number of stable vacua. The domain 
wall separating the two true vacua can be regarded in this case as the fixed plane 
under Z 2 transformations of the coordinate transverse to the plane. The three-form 
field Cn Ua will be able to couple to this wall in a manner discussed in the previous 
subsection. 



6 Axions and vacuum structure in gluodynamics 

6.1 Two scenarios 

The invisible axion is very light. Integrating out all other degrees of freedom and 
studying the low-energy axion effective Lagrangian must be a good approximation. 
The axion effective potential in QCD can be of two distinct types. 

Assuming that for all values of 9 the QCD vacuum is unique one arrives at the 
axion effective Lagrangian of the form 

-(d^a) 2 + ml (cos(a + 9) — 1] 



C a = f 



(79) 



2 

The axion potential does not have to be (and generically is not) a pure cosine; it 
may have higher harmonics. In the general case it is a smooth periodic function 
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of a + 9, with the period 2ir. For illustration we presented the potential as a pure 
cosine. This does not change the overall picture in the qualitative aspect. 

As we will see below, a smooth effective potential of the type ( [79"D emerges even if 
the (hadronic) vacuum family is non-trivial, but the transition between the distinct 
hadronic vacua does not occur inside the axion wall. This is the case with very light 
quarks, m q <C A/N c . In the opposite limit, one arrives at the axion potential with 
cusps, considered below. 

In the theory ([79]) one finds the axion walls interpolating between the vacuum 
state at a = — 9 and the same vacuum state at a = — 9 + 2n, 

a(z) + # = 4arctan(e ma2 ), (80) 

where the wall is assumed to lie in the xy plane, so that the wall profile depends 
only on z. This is the most primitive "27r wall." 

The tension of this wall is obviously of the order of 

Tx ~ f 2 m a . (81) 

Taking into account that f 2 m 2 ~ X where X is the topological susceptibility of the 
QCD vacuum, we get 

T x ~ X/m a . (82) 

The inverse proportionality to m a is due to the fact that the transverse size of the 
axion wall is very large. 

Let us now discuss the axion effective potential of the second type. In this case 
the potential has cusps, as is the case in pure gluodynamics, where the axion effective 
Lagrangian is of the form 

£a = y(^) 2 + min|^A 4 cos^±p} , (83) 

(see more detailed discussions below). Here the 9 angle was included in the definition 
of the axion field. The axion wall interpolates between a = and a = 2tt. 

What is the origin of this cusp? The cusps reflect a restructuring in the hadronic 
sector. When one (adiabatically) interpolates in a from to 2n a gluonic order 
parameter, for instance (GG), necessarily experiences a restructuring in the middle 
of the wall corresponding to the restructuring of heavy gluonic degrees of freedom. 
In other words, one jumps from the hadronic vacuum which initially (at a = 0) had 
(GG) = into the vacuum in which initially (GG) ^ 0. Upon arrival to a = 2tt, 
we find (GG) = again. This implies that the central part of such an axion wall 
is dominated by a gluonic wall. Thus, the cusp at a = 7r generically indicates the 



formation of a hadronic core, the D-wall |47j in the case at hand. 
Returning to the question of the tension we note that 

X ~ A 4 iy° , m a ~ in pure gluodynamics , 

X ~ A 3 N c m q , m a ~ A 3/2 ml /2 N^ 2 f' 1 in QCD with light quarks , (84) 
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which implies, in turn, 



fA 2 N® in pure gluodynamics 
T x ~ { (85) 
j A 3/2 m i/2 A r c 1 / 2 in QCD with light quarks . 

Here m q is the light quark mass. 

The presence of the large parameter / in T\ makes the axion halo the domi- 
nant contributor to the wall tension. The contribution of the hadronic component 
contains only hadronic parameters, although it may have a stronger dependence on 
N c . Examining the cusp with an appropriately high resolution one would observe 
that it is smoothed on the hadronic scale, where the hadronic component of the 
axion wall "sandwich" would become visible. The cusp carries a finite contribution 
to the wall tension which cannot be calculated in the low-energy approximation 
89]. To this end one needs to consider the hadronic core explicitly. The tension 



of the core T core ~ A 3 iV c , while the tension of the axion halo Th a i Q ~ fA 2 (in pure 
gluodynamics). 

We pause here to make a comment on the literature. The consideration of the 
axion walls in conjunction with hadrons dates back to the work of Huang and Sikivie, 
see Ref. ||46|| . This work treats the Weinberg- Wilczek N = 2 axion in QCD with 
two light flavors, which is replaced by a chiral Lagrangian for the pions, to the 
leading order (quadratic in derivatives and linear in the light quark masses). It is 
well-known [[361 , [371 |4lf that in this theory the crossover phenomenon takes place at 
m u = ma- In the realistic situation, (ma — m u ) / (ma + m u ) ~ 0.3 considered in Ref. 
||46|| . there is no crossover. The pions can be integrated over, leaving one with an 
effective Lagrangian for the axion of the type ( |75| ) (with a — > 2a). The potential 
is not pure cosine, higher harmonics occur too. The axion halo exhausts the wall, 
there is no hadronic core in this case. 

At the same time, Huang and Sikivie (see Ref. |46j) found an explicit solution for 



the "7T°" component of the wall. In fact, this is an illusion. The Huang-Sikivie (HS) 
solution refers to the bare ir° field. To find the physical 7r° field one must diagonalize 
the mass matrix at every given value of a (the bare fa is the physical axion field 
up to small corrections ~ f 2 /f 2 where /„- stands for the pion decay constant). Once 
this is done, one observes that the physical pion field, which is a combination of the 
bare pion and fa, is not excited in the HS solution. The equation (2.16) in the HS 
paper is exactly the condition of vanishing of the physical pion in the wall profile. 
This explains why the wall thickness in the HS work is of order m~ l , with no traces 
of the m" 1 component. The crossover of the hadronic vacua at a = ir/2 (remember, 
this is N = 2 model) could be recovered in the Huang-Sikivie analysis at m u = ma- 
However, the chiral pion Lagrangian predicts in this case the vanishing of the pion 
mass in the middle of the wall, for accidental reasons. This is explained in detail by 
A. Smilga, Ref. [M . 
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6.2 An Illustrative model 



To find the axion walls with D-wall core one has to solve QCD, which is way beyond 
our possibilities. Our task is more modest. We would like to obtain a qualitative de- 
scription of the axion wall sandwich which, with luck, can become semi- quantitative. 
To this end we want to develop toy models. An obvious requirement to any toy model 
is that it must qualitatively reproduce the basic features of the vacuum structure 
which we expect in QCD. In SUSY gluodynamics it was possible to write down a 
toy model with a Zjy c symmetry |8TJ which "integrates in" the heavy degrees of 
freedom and allows one to investigate the BPS domain walls in the large N c limit 
f54fl (see also [p8|). We will suggest a similar model in (nonsupersymmetric) QCD, 
then switch on axions, and study the axion domain walls in a semi-realistic setting. 
In this model we will be able to find exact solutions for D-walls and axion walls. 

Here we suggest a simple toy model which has a proper vacuum structure. In 
what follows an appropriate (complex) glue order parameter is denoted by $. The 
modulus and phase of this field describe respectively the ++ and 0~ + channels of 
the theory. 

Our toy model Lagrangian is 



C 

Vo 

Vi 



N 2 C A 2 



1 + 



Nr. 



'1 - $ N °e 



V = V + V 1 



+ 



XNl 



+ H.c. 



36) 



Here A is a numerical constant of order one, and X is the vacuum topological 
susceptibility in pure gluodynamics (note that X is independent of N c ). The scale 
parameter A is set to unity. 

This model has the vacuum family composed of N c states. Indeed, the minima 
of the energy are determined from the equations 



dV 



dV 







37) 



which have the following solutions (we put temporarily A = 1): 

9 + 2iri\ 



^vac 



exp i 



N c 



0,l,...,iV c -l 



(88) 



In the £-th minimum Vo vanishes, while VI produces a non- vanishing vacuum energy 
density, 



E t = XN 2 \l - cos 



9 + 2tt£ 

N r . 



(89) 
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For each given 9 the genuine vacuum is found by minimization, 



E{9) = N'Xmm, jl - cos (?±^1) ) . (90) 



* c 



The remaining N c —1 minima are quasivacua. Once the heavy field $ is integrated 
out, the vacuum energy is given by the expression (p0|); it has cusps at 9 = tt,3it 
and so on. Needless to say that the potential (|86|) has no cusps. 

We will first consider the model flSEp without the axion field, at 9 = 0, in the 
limit N c — > oo. In this limit the false vacua from the vacuum family are stable. 

The classical equation of motion defining the wall is 

^*" = §, (91) 

where primes denote differentiation with respect to z (we look for a solution which 
depends on the z coordinate only). 

This is a differential equation of the second order. It is possible, however, to 
reduce it to a first order equation. Indeed, Eq. (|9T| ) has an obvious "integral of 
motion" ("energy"), 

N 2 C $*'$' - V = Const = , (92) 

where the second equality follows from the boundary conditions. In the large N c 
limit one can parametrize the field $ as follows (p ~ 1): 

+ (93) 



Taking the square root of Eq. (p2|), substituting Eq. (|9~3"D and neglecting the terms 
of the subleading order in 1/N C we arrive at 

p = iAN c (1 - expp) . (94) 



The phase on the right-hand side can be chosen arbitrarily. The choice in Eq. (|94|) 
is made in such a way as to make it compatible with the boundary conditions for 
the wall interpolating between $ vac = 1 and $ vac = exp(27ri/A^ c ). This is precisely 
the expression that defines the domain walls in SUSY gluodynamics [0, |B8[] . It is 
not surprising that the same equation determines the D-walls in non-SUSY gluody- 
namics - the fermion-induced effects are not important for D-walls in the large N c 
limit. 



The solution of this equation was obtained in ||88|| . In the parametrization p 
a + ir the solution takes the form: 



cost = (p + 1) exp(— a 

ra{z) 
r(0) 



r (2) [exp(2t) - (1 + t) 2 }~ l/2 dt = -AN c \z\ . (95) 

J<t(0) 
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The real part of p is a bell-shaped function with an extremum at zero; it vanishes 
at ±00. The imaginary part of p, on the other hand, changes its value from to 
2tt. This determines a D-wall in the large N c gluodynamics. The width of the wall 
scales as l/N c . 

The solution presented above is exactly the same as in SUSY gluodynamics. 
This is not surprising since the ansatz fl93|) implies that V\ does not affect the 
solution - its impact is subleading in 1/N C , while Vq is exactly the same as in the 
SUSY-gluodynamics-inspired model of Ref. ||54j| . Moreover, for the same reason the 
domain wall junctions emerging in this model will be exactly the same as in the 
SUSY-gluodynamics-inspired model ||47||. Inclusion of V\ in the subleading order 



makes the wall to decay. 

Inclusion of the N = 1 axion field amounts to the replacement 

9 -> 6 + a 

in Eq. flSB]), plus the axion kinetic term 

£ kin = ( - + 2 2$ *^ (d,af + iN c (d„a)($*d^ - $.9^*) . (96) 

The occurrence of the mixing between a and the phase of $ is necessary, as is readily 
seen from the softly broken SUSY gluodynamics. (To get the potential of the type 
(|56j) in this model, one must eliminate the GG term by a chiral rotation. Then 
m — > m exp((# + a) / N c ) and, additionally one gets d^ax [the gluino axial current].) 
The term 2$*$ in the brackets has to be included to reproduce the correct mass 
for the axion after the physical heavy state is integrated out. The presence of this 
term signals that QCD dynamics generates not only the potential for the axion but 
also modifies its kinetic term. On the other hand, since $*$ < A 2 and, moreover, 
A << /, this term can be neglected for all practical purposes. 

We are interested in the configuration with a interpolating between and 2n. 
The phase of $ will first adiabatically follow a/N c , then at a ~ it, when the phase 
of $ is close to tt/N c , it will very quickly jump by —27r/N c , and then it will continue 
to grow as a/N c , so that when a reaches 2ir the phase of $ returns to zero. This 
jump is continuous, although it occurs at a scale much shorter than m~ x . This 
imitates the D-wall core of the axion wall. One cannot avoid forming this core, 
since otherwise the interpolation would not connect degenerate states - on one side 
of the wall we would have (hadronic) vacuum, on the other side an excited state. 

In the large N c limit one can be somewhat more quantitative. Indeed, in this 
approximation the model admits the exact solutions. The gluonic core of the wall 
has the same form as before, Eq. (0), but the phase r is now substituted by the 
superposition r — (a + 6) since the axion field is mixed with the phase of the <3? field. 

This very narrow core is surrounded by a diffused axion halo. The axion field is 
described in this halo by the solution to the Lagrangian (^). This takes the form: 

6 + a(z) = -27T + 4N C arctan (e maZ tan-^-j , z < 0, 
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e + a(z) = - AN C arctan [e~ maZ tan— J , z > . (97) 

Thus, we find explicitly the stable axion wall with a D-wall core. Note that this is 
a usual "27r" wall as it separates two identical hadronic vacua. As we discussed in 
the introduction, this wall can decay quantum mechanically. However, its lifetime 
is infinite for all practical purposes. Moreover, this wall is harmless cosmologically. 
It will be produced bounded by global axion strings in the early universe. Bounded 
walls shrink very quickly by decaying into axions and hadrons. 



7 Axions and vacuum structure in QCD with light 
quarks 

So far we discussed pure gluodynamics with the axion. Our final goal is to study 
QCD with Nf = 3. There are two, physically distinct regimes to be considered in 
this case. In real QCD 

m u , m d < m s ~ , m u , m d , m s < A . (98) 

In this regime the consideration of the chiral Lagrangians |3^, |3?|, |4"I|| , does not 
exhibit the vacuum family. We will comment on why the light quarks screen the 
vacuum family of the glue sector, so that the axion domain wall provides no access 
to it. In the limit (p8|) the effects due to the D-walls will be marginal. 
On the other hand, in the genuinely large N c limit 

A . . , 

— < m u , m d < m s < A, (99) 

physics is rather similar to that of pure gluodynamics. The light quarks are too 
heavy to screen the vacuum family of the glue sector 

In what follows we study the axion walls and their hadronic components in the 



limits (|98|) and (0), separately. 



7.1 One light quark 

To warm up, let us start from the theory with one light quark. In the limit of large 
A^ c this introduces a light meson, "r/"' . An appropriate effective Lagrangian can 
be obtained by combining the vacuum energy density of gluodynamics with what 
remains from the chiral Lagrangian at Nf = 1, 



2 



V = -m q K A N c cos /3 + min <j -A^ C 2 A 4 cos - + 9 J~ ^ \ . (100) 



N c 
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Here (3 is the phase of U ~ QlQr, while F 2 ~ A 2 A^ C is the "t/" coupling constant 
squared. The product Ff3 is the "77"' field. The first term in V corresponds to the 
quark mass term, A4U + h.c. At N c = 00 the second term in V becomes ((3 + 9) 2 . It 
corresponds to (iln detU + 9) 2 in Eq. (11) in [36]. The subleading in 1/N C terms sum 
up into a 2n periodic function of the cosine type, with the cusps. It is unimportant 
that we used cosine in Eq. ( |100| ). Any 2ir periodic function of this type would lead 



to the same conclusions. The second term in Eq. (|100|) differs from the vacuum 
energy density in gluodynamics by the replacement 9 — > (3 + 9. 

If m q <C A/N c , the first term in V is a small perturbation; therefore, in the 
vacuum, (3 + 9 = 2irk, and, hence, the 9 dependence of the vacuum energy is 

£ vac (9) = -m q A 3 N c cos9. (101) 

It is smooth, 2ir periodic and proportional to m q as it should be on general grounds 
in the theory with one light quark. 

The condition m q C A/Af c precludes us from sending N c — > 00. The would be 
"27r" wall in the variable (3 is expected to be unstable. This is due to the fact that 
at N c ~ 3 the absolute value of the quark condensate ipip is not "harder" than the 
phase of the condensate (3, and the barrier preventing the creation of holes in the 
"27r" wall is practically absent. 

If one closes one's eyes on this instability one can estimate that the tension of 
the "77"' wall is proportional to A 3 N^ 2 , with a small correction m q A 2 N 3 ^ 2 from the 
quark mass term. The tension of the D wall core is, as previously, A 3 N C . 

In the opposite limit 



m q > — , but m q still < A , (102) 



the situation is trickier. Now the first term in V is dominant, while the second is a 
small perturbation. There are A^ c distinct vacua in the theory, 

2A 

= (0 + 2-ni) . (103) 

m q N c 

Then the 9 dependence of the vacuum energy density is 

S vac (9) = A 4 min(0 + 2tv£) 2 , (104) 

this is similar to that in the theory without light quarks (i.e., the same as in gluo- 
dynamics). The "77"' wall is stable at A^ c — > 00, with a a D-wall core in its center. 
The rj' wall is a u 2ir" wall. 

From this standpoint the quark with the mass ( |102j ) is already heavy, although 
the "77"' is still light on the scale of A, 

M^~mfA 1 / 2 «A. 
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So far the axion was switched off. What changes if one includes it in the theory? 
The Lagrangian now becomes 



£ = ^(^) 2 + ^(«vo 2 -m«), 

V = —m q A 3 N c cos (3 + min j - iV c 2 A 4 cos ^ + ^ + ^ j , (105) 

where the 9 angle is absorbed in the definition of the axion field. 

The bare "r/" mixes with the bare axion. It is easy to see that in the limit 
m q <C A/N c the physical "77"' is proportional to (3 + a, rather than to /3. Therefore, 
even if we force the axion wall to develop, (i.e. a to evolve from to 2ir) the "77"' 
wall need not develop. It is energetically expedient to have (3 + a = 0. Thus, the 
effect of the axion field on the hadronic sector is totally screened by a dynamical 
phase f3 coming from the quark condensate. In other words, the axion wall with the 
lowest tension corresponds to the frozen physical "77"', 

(3 + a = 0. 

There is no hadronic core. The tension of this wall is determined from the term 
oc m q A 3 N c . 

(If one wishes, one could add the "77"' wall to the axion wall. Then the "77"' 
wall, with the D-wall core will appear in the middle of the axion wall, but they are 
basically unrelated. This will be a secondary phenomenon, and the D wall core will 
be, in fact, the core of the "77"' wall rather than the axion wall.) 

If the quark mass is such that ( 102 ) applies, then the axion field a cannot be 



screened, since we cannot freeze (3-\-a everywhere in the axion wall profile at zero - 
at m q ^> A/N c , (3 is proportional to the physical "77"' and is much heavier than the 
axion field. Thus, in this case the axion wall will be described by the Lagrangian 
( p3l) and will have a D-wall core. One may also add, on top of it, the "77"' wall. This 
will cost m*' 2 A 5 ' 2 Af c in the wall tension - still much less than A 3 iV c of the D-wall 
core of the axion wall. 

The limit ( |102| ) is unrealistic. Moreover, in this limit the D walls taken in iso- 
lation, without the axion walls, are stable by themselves, although they interpolate 



between nondegenerate states J45]. 



7.2 Three Light Quarks 

Let us turn the case of three light flavors. The physical picture is quite similar to 
that of the one-flavor case, see Sec. 7.1. 

We assume the mass matrix M. in the meson Lagrangian to be diagonal. There- 
fore, we will looking for a diagonal U(3) meson matrix which minimizes the potential, 

U = diag (e^ 1 , e^ 2 , e i03 ) . (106) 
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The potential takes the form 



V = -J^m^N.cos^ + min \-N 2 c A 4 cos gi^i±^±j^ | . (107 ) 



* c 



As before, we will consider two limiting cases, ( |98| ) and (p9|). 

Let us switch off the axion field first. In the limit of genuinely light quarks, Eq. 
(PSf), when the second term in the potential (|107|) is dominant, the solutions for 0's 



were found in [56 , 57| . They satisfy to the relation 3 ~ and 0i + 02 = — The 
corresponding expression for the vacuum energy density is 



£vac(#) = -N c A d ^/ ml + m 2 d + 2m u m d cos 9 . (108) 

As in Sec. 7.1, we deal here with a smooth single- valued function of 9. The inclusion 
of the axion replaces 9 — > 9 + a — > a. The physical rj' field is given by the sum 
J2i (pi + a - It is energetically favorable to freeze this state. Thus, the situation is 
identical to that in the one-flavor case: even if the axion wall is forced to develop, the 
physical rj' wall (which is now the J2i 4>i + a wall) does not have to occur. Effectively, 
the vacuum angle is screened, and there is no D-wall core in the axion wall. 

If, nonetheless, the rj' wall is formed due to some cosmological initial conditions, it 
will have a D-wall core (albeit the rj' wall is unstable in the limit at hand and cannot 
be considered in the static approximation). The would-be rj' wall is independent of 
the axion wall; its effect on the axion wall formation is rather irrelevant. 

In addition to this, a "27r" wall could develop for nonsinglet mesons at certain 
values of the quark masses. There is nothing new we could add to this issue which 
is decoupled from the issue of the vacuum family in the glue sector and D-walls. 

We now pass to the opposite limit fl9^), when the first term in the potential (|107| ) 
is dominant. As in Sec. 7.1, there are N c distinct vacua with the energy given by 
( |104|) . It is straightforward to show that the potential for the axion in this case is 
of the form (|53"D, with the cusps which signal the presence of the D-wall core. This 
is similar to what happens in gluodynamics. One cannot avoid having an 7/ wall 
in the middle of the axion wall, which entails a D-wall too. The D-walls separate 
the degenerate vacua. Since they "live" in the middle of the axion wall, they are 
perfectly stable. 

(In addition, there can be "27r" walls in either of 0's or their linear combinations. 
However, these latter are unstable and do not appear in the physical spectrum of 
the theory.) 
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